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Abstract.  We  compare  Stetson’s  photometric  standards  with  measurements 
listed  in  a  standard  star  catalog  constructed  using  repeated  SDSS  imaging  ob¬ 
servations.  The  SDSS  catalog  includes  over  700,000  candidate  standard  stars 
from  the  equatorial  stripe  82  ( |  DEC  |  <  1?266  )  in  the  RA  range  20h  34m  to  4h 
00m,  and  with  the  r  band  magnitudes  in  the  range  14-21.  The  distributions  of 
measurements  for  individual  sources  demonstrate  that  the  SDSS  photometric 
pipeline  correctly  estimates  random  photometric  errors,  which  are  below  OUOl  for 
stars  brighter  than  (19.5,  20.5,  20.5,  20,  18.5)  in  ugriz,  respectively  (about  twice 
as  good  as  for  individual  SDSS  runs).  We  derive  mean  photometric  transfor¬ 
mations  between  the  SDSS  gri  and  the  BVRI  system  using  1165  Stetson  stars 
found  in  the  equatorial  stripe  82,  and  then  study  the  spatial  variation  of  the  dif¬ 
ference  in  zeropoints  between  the  two  catalogs.  Using  third-order  polynomials 
to  describe  the  color  terms,  we  find  that  photometric  measurements  for  main- 
sequence  stars  can  be  transformed  between  the  two  systems  with  systematic 
errors  smaller  than  a  few  millimagnitudes.  The  spatial  variation  of  photometric 
zeropoints  in  the  two  catalogs  typically  does  not  exceed  0!"  0 1 .  Consequently,  the 
SDSS  Standard  Star  Catalog  for  Stripe  82  can  be  used  to  calibrate  new  data  in 
both  the  SDSS  ugriz  and  the  BVRI  systems  with  a  similar  accuracy. 
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1.  Introduction 

Astronomical  photometric  data  are  usually  calibrated  using  sets  of  standard 
stars  whose  brightness  is  known  from  previous  work.  The  most  notable  modern 
optical  standard-star  catalogs  are  Landolt  standards  (Landolt  1992)  and  Stet¬ 
son  standards  (Stetson  2000,  2005).  Both  are  reported  on  the  Johnson-Kron- 
Cousins  system  (Landolt  1983  and  references  therein).  The  Landolt  catalog 
provides  magnitudes  accurate  to  1-2%  in  the  UBVRI  bands  for  ~500  stars  in 
the  V  magnitude  range  11.5-16.  Stetson  has  extended  Landolt’s  work  to  fainter 
magnitudes,  and  provided  the  community  with  ~l-2%  accurate  magnitudes  in 
the  BVRI  bands  for  ~15,000  stars  in  the  magnitude  range  V  20.  Most  stars 
from  both  sets  are  distributed  along  the  celestial  equator,  which  facilitates  their 
use  from  both  hemispheres. 

The  data  obtained  by  the  Sloan  Digital  Sky  Survey  (SDSS,  York  et  al. 
2000)  can  be  used  to  extend  the  work  by  Landolt  and  Stetson  to  even  fainter 
levels,  and  to  increase  the  number  of  standard  stars  to  almost  a  million.  In 
addition,  SDSS  has  designed  its  own  photometric  system  ( ugriz ,  Fukugita  et 
al.  1996)  which  is  now  in  use  at  a  large  number  of  observatories  worldwide.  This 
widespread  use  of  the  ugriz  photometric  system  motivates  the  construction  of  a 
large  standard  star  catalog  with  ~1%  accuracy.  As  a  part  of  its  imaging  survey, 
SDSS  has  obtained  many  scans  in  the  so-called  Stripe  82  region,  defined  by 
|DEC|  <  1?266  and  RA  approximately  in  the  range  20h-4h.  These  repeated 
observations  can  be  averaged  to  produce  more  accurate  photometry  than  the 
nominal  2%  single-scan  accuracy  (Ivezic  et  al.  2004). 

We  describe  the  construction  and  testing  of  a  standard  star  catalog  in  Sec.  2, 
and  discuss  its  comparison  with  Stetson’s  photometric  standards  in  Sec.  3. 


2.  The  SDSS  Standard  Star  Catalog  for  Stripe  82 
2.1.  Overview  of  SDSS  Imaging  Data 

SDSS  is  using  a  dedicated  2.5-m  telescope  (Gunn  et  al.  2006)  to  provide  homo¬ 
geneous  and  deep  (r  <  22.5)  photometry  in  five  passbands  (Fukugita  et  al.  1996; 
Gunn  et  al.  1998;  Smith  et  al.  2002;  Hogg  et  al.  2002)  repeatable  to  01n02  (rms 
scatter,  hereafter  rms,  for  sources  not  limited  by  photon  statistics,  Ivezic  et 
al.  2003)  and  with  a  zeropoint  uncertainty  of  ~0. 02-0. 03  (Ivezic  et  al.  2004). 
The  survey  sky  coverage  of  close  to  ~10,000  deg2  in  the  Northern  Galactic  Cap, 
and  ~300  deg2  in  the  southern  Galactic  hemisphere,  will  result  in  photometric 
measurements  for  well  over  100  million  stars  and  a  similar  number  of  galaxies1. 
Astrometric  positions  are  accurate  to  better  than  O'!  1  per  coordinate  (rms)  for 
sources  with  r  <  20“15  (Pier  et  al.  2003),  and  the  morphological  information 
from  the  images  allows  reliable  star-galaxy  separation  to  r  ~  21.5  (Lupton  et 
al.  2001,  2003;  Scranton  et  al.  2002). 

Data  from  the  imaging  camera  (thirty  photometric,  twelve  astrometric,  and 
two  focus  CCDs,  Gunn  et  al.  1998)  are  collected  in  drift  scan  mode.  The 


lrThe  recent  Data  Release  5  lists  photometric  data  for  215  million  unique  objects  observed  in 
8000  deg2  of  sky,  please  see  http://www.sdss.org/dr5/ 
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images  that  correspond  to  the  same  sky  location  in  each  of  the  five  photometric 
bandpasses  (these  five  images  are  collected  over  ~5  minutes,  with  54  sec  for  each 
exposure)  are  grouped  together  for  processing  as  a  field.  A  field  is  defined  as 
a  36  seconds  (1361  pixels,  or  9  arcmin)  long  and  2048  pixels  wide  (13  arcmin) 
stretch  of  drift-scanning  data  from  a  single  column  of  CCDs  (sometimes  called 
a  scanline,  for  more  details  please  see  Stoughton  et  al.  2002,  Abazajian  et  al. 
2003,  2004,  2005,  Adelman- McCarthy  et  al.  2006).  Each  of  the  six  scanlines 
(called  together  a  strip)  is  13  arcmin  wide.  The  twelve  interleaved  scanlines  (or 
two  strips)  are  called  a  stripe  (~2?5  wide). 

2.2.  The  Photometric  Calibration  of  SDSS  Imaging  Data 

SDSS  2.5-m  imaging  data  are  photometrically  calibrated  using  a  network  of 
calibration  stars  obtained  in  1520  41.5  x  41.5  arcmin2  transfer  fields,  called  sec¬ 
ondary  patches.  These  patches  are  positioned  throughout  the  survey  area  and 
are  calibrated  using  a  primary  standard-star  network  of  158  stars  distributed 
around  the  northern  sky  (Smith  et  al.  2002).  The  primary  standard  star  network 
is  tied  to  an  absolute  flux  system  by  the  single  F0  subdwarf  star  BD+174708, 
whose  absolute  fluxes  in  SDSS  bands  are  taken  from  Fukugita  et  al.  (1996). 
The  secondary  patches  are  grouped  into  sets  of  four,  and  are  observed  by  the 
Photometric  Telescope  (hereafter  PT)  in  parallel  with  observations  of  the  pri¬ 
mary  standards.  A  set  of  four  spans  all  12  scanlines  of  a  survey  stripe  along 
the  width  of  the  stripe,  and  the  sets  are  spaced  along  the  length  of  a  stripe  at 
roughly  15°  intervals,  which  corresponds  to  an  hour  of  scanning  at  the  sidereal 
rate. 

SDSS  2.5-nr  magnitudes  are  reported  on  the  “natural  system”  of  the  2.5-m 
telescope  defined  by  the  photon-weighted  effective  wavelengths  of  each  combi¬ 
nation  of  SDSS  filter,  CCD  response,  telescope  transmission,  and  atmospheric 
transmission  at  a  reference  airmass  of  1.3  as  measured  at  APO.2  The  magnitudes 
are  referred  to  as  the  ugriz  system  (which  is  different  from  the  “primed”  system, 
u'g'r'i'z',  that  is  defined  by  the  PT3).  The  reported  magnitudes4  are  corrected 
for  atmospheric  extinction  (using  simultaneous  observations  of  standard  stars 
by  the  PT)  and  thus  correspond  to  measurements  at  the  top  of  the  atmosphere'5 
(except  for  the  fact  that  the  atmosphere  has  an  impact  on  the  wavelength  de¬ 
pendence  of  the  photometric  system  response).  The  magnitudes  are  reported  on 
the  AB  system  (Oke  &;  Gunn  1983)  defined  such  that  an  object  with  a  specific 
flux  of  Fj,=3631  Jy  has  m  =  0  (i.e.  an  object  with  Fj,=constant  has  an  AB 
magnitude  equal  to  the  Johnson  V  magnitude  at  all  wavelengths).  In  summary, 


2  Transmission  curves  for  the  SDSS  2.5-m  photometric  system  are  available  at 
http : //www. sdss . org/dr5/ instruments/ imager/ 

3For  subtle  effects  that  led  to  this  distinction,  please  see  Stoughton  et  al.  (2002)  and 
http : //www. sdss . org/dr5/algorithms/f luxcal .html 

4SDSS  uses  a  modified  magnitude  system  (Lupton  et  al.  1999),  which  is  virtually  identical  to 
the  standard  astronomical  Pogson  magnitude  system  at  high  S/N  ratios  relevant  here. 

®The  same  atmospheric  extinction  correction  is  applied  irrespective  of  the  source  color;  the 
systematic  errors  this  introduces  are  probably  less  than  1%  for  all  but  objects  of  the  most 
extreme  colors. 
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given  a  specific  flux  of  an  object  at  the  top  of  the  atmosphere,  Fv( A),  the  re¬ 
ported  SDSS  2.5-rn  magnitude  in  a  given  band  corresponds  to  (modulo  random 
and  systematic  errors,  which  will  be  discussed  later) 

m  =  -2.5i°gltl(gA^),  (1) 

where 

F0  =  J  Fu(X)(j)(X)dX.  (2) 

Here,  0(A)  is  the  normalized  system  response  for  the  given  band, 


0(A) 


A_1<S(A) 

fX-1S(X)dX’ 


(3) 


with  the  overall  atmosphere+system  throughput,  S( A),  available  from  the  web¬ 
site  given  above  (see  Smolcic  et  al.  2006  for  a  figure  showing  0(A)  for  the  ugriz 
system) . 

The  quality  of  SDSS  photometry  stands  out  among  available  large-area 
optical  sky  surveys  (Ivezic  et  al.  2003,  2004;  Sesar  et  al.  2006).  Nevertheless, 
the  achieved  accuracy  is  occasionally  worse  than  the  nominal  0“02-0™03  (rms 
scatter  for  sources  not  limited  by  photon  statistics).  Typical  causes  of  substan¬ 
dard  photometry  include  an  incorrectly  modeled  PSF  (usually  due  to  fast  varia¬ 
tions  of  atmospheric  seeing,  or  lack  of  a  sufficient  number  of  the  isolated  bright 
stars  needed  for  modeling  the  PSF),  unrecognized  changes  in  atmospheric  trans¬ 
parency,  errors  in  photometric  zeropoint  calibration,  effects  of  crowded  fields 
at  low  galactic  latitudes,  an  undersampled  PSF  in  excellent  seeing  conditions 
(<;  0'/8;  the  pixel  size  is  (f/4),  incorrect  flatfield,  or  bias  vectors,  scattered  light 
correction,  etc.  Such  effects  can  conspire  to  increase  the  photometric  errors  to 
levels  as  high  as  0H05  (with  a  frequency,  at  that  error  level,  of  roughly  one  field 
per  thousand).  However,  when  multiple  scans  of  the  same  sky  region  are  avail¬ 
able,  many  of  these  errors  can  be  minimized  by  properly  averaging  photometric 
measurements. 


2.3.  The  Catalog  Construction  and  Internal  Tests 

A  detailed  description  of  the  catalog  construction,  including  the  flatfield  correc¬ 
tions,  and  various  tests  of  its  photometric  quality  can  be  found  in  Ivezic  et  al. 
(2006).  Here  we  only  briefly  describe  the  main  catalog  properties. 

The  catalog  is  based  on  58  SDSS  runs  from  stripe  82  (approximately  20h 
<  RA  <  04h  and  |DEC|  <  1?266)  obtained  in  mostly  photometric  conditions 
(as  indicated  by  the  calibration  residuals,  infrared  cloud  camera6,  and  tests 
performed  by  the  runQA  quality  assessment  pipeline7).  Candidate  standard  stars 
from  each  run  are  selected  by  requiring 

1.  that  objects  are  unresolved  (classified  as  STAR  by  the  photometric  pipeline), 


°For  more  details  about  the  camera  see  http://hoggpt.apo.nmsu.edu/irsc/irsc_doc/ 
7For  a  description  of  runQA  pipeline,  see  Ivezic  et  al.  2004. 
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2.  that  they  have  quoted  photometric  errors  (as  computed  by  the  photometric 
pipeline)  smaller  than  0™05  in  at  least  one  band,  and 

3.  that  processing  flags  BRIGHT,  SATUR,  BLENDED,  EDGE  are  not  set8. 

These  requirements  select  unsaturated  sources  with  sufficiently  high  signal-to- 
noise  per  single  observation  to  approach  the  final  photometric  errors  of  0™02  or 
smaller. 

After  positionally  matching  (within  1  arcsec)  all  detections  of  a  single  star, 
various  photometric  statistics  such  as  mean,  median,  rms  scatter,  number  of 
observations,  and  y2  per  degree  of  freedom  (y2df)  are  computed  in  each  band. 
This  initial  catalog  of  multi-epoch  observations  includes  924,266  stars  with  at 
least  4  observations  in  each  of  the  g ,  r  and  i  bands.  The  median  number  of 
observations  per  star  and  band  is  10,  and  the  total  number  of  photometric 
measurements  is  ~40  million.  The  errors  for  the  averaged  photometry  are  below 
CP01  at  the  bright  end.  These  errors  are  reliably  computed  by  the  photometric 
pipeline,  as  indicated  by  the  y2df  distributions. 

Adopted  candidate  standard  stars  must  have  at  least  4  observations  in  each 
of  the  g,  r  and  i  bands  and,  to  avoid  variable  sources,  y2df  less  than  3  in  the  gri 
bands.  The  latter  cut  rejects  about  20%  of  stars.  We  also  limit  the  RA  range  to 
20h  34m  <  RA  <  04h  00m,  which  provides  a  simple  areal  definition  (together 
with  |DEC|  <1?266)  while  excluding  only  a  negligible  fraction  of  stars.  These 
requirements  result  in  a  catalog  with  681,262  candidate  standard  stars.  Of  those, 
638,671  have  the  random  error  for  the  median  magnitude  in  the  r  band  smaller 
than  0r.n01,  and  131,014  stars  have  the  random  error  for  the  median  magnitude 
smaller  than  CP01  in  all  five  bands. 

The  internal  photometric  consistency  of  this  catalog  is  tested  using  a  variety 
of  methods,  including  the  position  of  the  stellar  locus  in  the  multi-dimensional 
color  space,  color-redshift  relations  for  luminous  red  galaxies,  and  a  direct  com¬ 
parison  with  the  secondary  standard  star  network.  While  none  of  these  methods 
is  without  its  disadvantages,  together  they  suggest  that  the  internal  photometric 
zeropoints  are  spatially  stable  at  the  1%  level. 

3.  An  External  Test  of  Catalog  Quality  Based  on  Stetson’s  Standards 

While  the  tests  based  on  SDSS  data  suggest  that  the  internal  photometric  ze¬ 
ropoints  are  spatially  stable  at  the  1%  level,  it  is  of  course  prudent  to  verify 
this  conclusion  using  an  external  independent  dataset.  The  only  large  external 
dataset  with  sufficient  overlap,  depth  and  accuracy  to  test  the  quality  of  the 
Stripe  82  catalog  is  that  provided  by  Stetson  (2000,  2005).  Stetson’s  catalog 
lists  photometry  in  the  BVRI  bands  (Stetson’s  photometry  is  tied  to  Landolt’s 
standards)  for  ~1,200  stars  in  common  (most  have  V  <  19.5).  We  synthesize 
the  BVRI  photometry  from  SDSS  gri  measurements  using  photometric  trans¬ 
formations  of  the  following  form 

^Stetson  -  MSDSS  =  A  C3  +  B  C2  +  C  C  +  D,  (4) 


8  For  more  details  about  photometric  processing  flags  see  Stoughton  et  al.  (2002)  and 
http : //www. sdss . org/dr4/products/ catalogs/f lags .html 
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Figure  1.  An  illustration  of  the  need  for  non-linear  color  terms  when  trans¬ 
forming  SDSS  photometry  to  the  BVRI  system.  The  small  dots  in  the  top 
panel  show  the  residuals  (in  magnitudes)  for  the  cubic  B  —  g  transformation 
based  on  Eq.  4,  as  a  function  of  the  g—r  color.  The  large  symbols  show  the  me¬ 
dians  for  Ob1!  wide  g  —  r  bins.  The  second  panel  is  analogous  to  the  top  panel, 
except  that  a  best-fit  linear  transformation  is  used  B—g  =  0.345  (g— r)+0.205. 
Note  the  increased  deviation  of  median  residuals  from  zero.  The  third  panel 
shows  residuals  for  the  relation  B  —  V  =  0.949  (g  —  r)  +  0.197  and  demon¬ 
strates  that  even  such  a  color  vs.  color  relation  is  measurably  non-linear.  The 
bottom  panel  shows  that  this  is  not  a  peculiarity  of  the  g—  color  because  a 
transformation  based  on  the  B  —  V  color,  B  —  g  =  0.364  (B  —  V)  +  0.133,  is 
also  measurably  non-linear. 


where  m  =  (BVRI)  and  g  =  (g,  g ,  r,  z),  respectively,  and  the  color  c  is  measured 
by  SDSS  (g  —  r  for  the  B  and  V  transformations,  and  r  —  i  for  the  R  and  / 
transformations).  The  measurements  are  not  corrected  for  the  ISM  reddening. 
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Table  1.  SDSS  to  BVRI  transformations. 


color 

med 

^med 

Xmed 

<  >  all 

^all 

A 

B 

C 

D 

B-g 

-1.6 

8.7 

1.4 

1.0 

32 

0.2628 

-0.7952 

1.0544 

0.0268 

V-g 

0.8 

3.9 

1.0 

0.9 

18 

0.0688 

-0.2056 

-0.3838 

-0.0534 

R—  r 

-0.1 

5.8 

0.9 

1.2 

15 

-0.0107 

0.0050 

-0.2689 

-0.1540 

I  -  i 

0.9 

6.1 

1.0 

1.2 

19 

-0.0307 

0.1163 

-0.3341 

-0.3584 

<>med •  the  median  value  of  median  transformation  residuals  (differences  between  the  measured  values 
of  colors  listed  in  the  first  column  and  those  synthesized  using  Eq.  4)  in  0^1  wide  g  —  r  bins  for  stars  with 
0.25<  g  —  r  <1.45  (in  millimag).  These  medians  of  medians  measure  the  typical  level  of  systematics  in 
the  gri-to-BVRI  photometric  transformations  introduced  by  the  adopted  analytic  form  (see  Eq.  4). 
(7raed:  the  rms  scatter  for  median  residuals  described  above  (in  millimag). 

Xmed  •  the  rms  scatter  for  residuals  normalized  by  statistical  noise.  The  listed  values  are  ~1,  which 
indicates  that  the  scatter  around  adopted  photometric  transformations  listed  under  crme(j  is  consistent 
with  expected  noise. 

<>aii:  the  median  value  of  residuals  evaluated  for  all  stars  (in  millimag). 
crall :  the  median  value  of  residuals  evaluated  for  all  stars  (in  millimag). 
e:  the  rms  scatter  for  residuals  evaluated  for  all  stars  (in  millimag). 

f:  coefficients  A-D  needed  to  transform  SDSS  photometry  to  the  BVRI  system  (see  Eq.  4). 


Traditionally,  such  transformations  are  assumed  to  be  linear  in  color9.  We  use 
higher-order  terms  in  Eq.  4  because  at  the  1  2%  level  there  are  easily  detectable 
deviations  from  linearity  for  all  color  choices,  as  shown  in  Fig.  1.  The  best- 
fit  coefficients  for  the  transformation  of  SDSS  gri  measurements  to  the  BVRI 
system10  are  listed  in  Table  1,  as  well  as  low-order  statistics  for  the  mgtetson  — 
/iSDSS  difference  distribution.  We  find  no  trends  as  a  function  of  magnitude  at 
the  <  0E005  level. 

We  have  also  tested  for  the  effects  of  interstellar  dust  reddening  and  rnetal- 
licity  on  the  adopted  photometric  relations.  For  about  half  of  stars  in  common, 
the  SFD  map  (Schlegel,  Finkbeiner  Sz  Davis  1998)  lists  Eb~v  >  0.15.  The  dif¬ 
ferences  in  median  residuals  for  these  stars  and  those  with  smaller  Eb~v  (the 
median  Eb~v  are  0.31  and  0.04)  are  always  less  than  0™01  (the  largest  difference 
is  0™008  for  the  B  —  g  transformation). 

Stars  at  the  blue  tip  of  the  stellar  locus  with  u  —  g  <  1  are  predominantly 
low-metallicity  stars  (Bond  et  al.  2006,  in  prep.),  as  illustrated  in  Fig.  2.  Fig.  3 
shows  that  the  median  residuals  for  mstetson  —  /isdss  are  the  same  for  the  0.8  < 
u  —  g  <  0.95  and  1.0  <  u  —  g  <  1.15  subsamples  to  within  their  measurement 
errors  (~0F010).  There  is  a  possibility  that  the  offset  is  somewhat  larger  for  the 
B  —  g  transformation  for  stars  with  u  —  g  <  0.9,  but  its  statistical  significance 
is  low.  If  this  is  a  true  effect,  it  implies  a  gradient  with  respect  to  metallicity  of 
about  0F02/dex. 

We  conclude  that  the  SDSS  catalog  described  here  could  also  be  used  to 
calibrate  the  data  to  the  BVRI  system  without  a  loss  of  accuracy  due  to  trans¬ 
formations  between  the  two  systems. 


9For  various  photometric  transformations  between  the  SDSS  and  other  systems,  see  Abazajian 
et  al.  (2005)  and  http://www.sdss.org/dr4/algorithms/sdssUBVRITransform.html 

10The  same  transformations  can  be  readily  used  to  transform  BVRI  measurements  to  the  corre¬ 
sponding  gri  values  because  B  —  V  =  f  (g  —  r)  and  R  —  I  =  f(r  —  i)  are  monotonous  functions. 
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median  metallicity  (-2.5  to  -0.5:  blue  to  red)  median  Temp.  (6700  K  to  5300  K:  blue  to  red) 


0.8  0.9  1  1.1  1.2  1.3  0.8  0.9  1  1.1  1.2  1.3 


u  —  g  u  —  g 

Figure  2.  The  median  effective  temperature  and  median  metallicity  esti¬ 
mated  from  SDSS  spectra  of  ~40,000  stars  by  Allende  Prieto  et  al.  (2006), 
shown  as  a  function  of  the  position  in  the  g  —  r  vs.  u  —  g  diagram  based  on 
SDSS  imaging  data.  In  the  left  panel,  the  temperature  in  each  color-color  bin 
is  linearly  color-coded  from  5300  K  (red)  to  6,700  K  (blue).  The  right  panel 
is  analogous,  except  that  it  shows  the  median  metallicity,  linearly  color-coded 
from  —0.5  (red)  to  —2.5  (blue).  Full-color  version  on  page  618. 


3.1.  The  Spatial  Variations  of  the  Zeropoints 


Table  2.  Photometric  zeropoint  spatial  variations. 


color 

XR1 

a Rl 

N 

ae 

to 

aR2 

N^2 

XR3 

aR3 

3 

W 

~.a 

XRA 

aR4. 

Nc 

R4 

B-g 

-29 

21 

92 

6 

27 

165 

8 

42 

155 

-4 

27 

281 

V-g 

0 

17 

99 

0 

15 

217 

6 

25 

161 

17 

19 

282 

R-r 

-6 

16 

58 

4 

16 

135 

-8 

12 

11 

39 

27 

60 

I  -  i 

-11 

16 

94 

6 

18 

205 

2 

16 

124 

19 

15 

47 

a:  The  median  value  of  residuals  (in  millimag)  for  transformations  listed  in  the  first  column,  evalu¬ 
ated  separately  for  regions  1-4,  defined  as:  Rl:  RA~325,  DEC<0;  R2:  RA~15;  R3:  RA~55;  R4: 
RA~325,  DEC>0. 

b:  The  rms  scatter  for  the  transformation  residuals  (in  millimagnitudes) . 

c:  The  number  of  stars  in  each  region  with  good  photometry  in  the  required  bands. 


The  BVRI  photometry  from  Stetson  and  that  synthesized  from  SDSS  agree 
at  the  level  of  0™02  (rms  scatter  for  the  magnitude  differences  of  individual 
stars;  note  that  the  systems  are  tied  to  each  other  to  within  a  few  millimags  by 
transformations  listed  in  Table  1).  This  scatter  is  consistent  with  the  claimed 
accuracy  of  both  catalogs  (the  magnitude  differences  normalized  by  the  implied 
error  bars  are  well  described  by  Gaussians  with  widths  in  the  range  0.7-0. 8). 
This  small  scatter  allows  us  to  test  for  the  spatial  variation  of  zeropoints  between 
the  two  datasets,  despite  the  relatively  small  number  of  stars  in  common. 

Stars  in  common  are  found  in  four  isolated  regions  that  coincide  with  histor¬ 
ical  and  well-known  Kapteyn  Selected  Areas  113,  92,  95,  and  113.  We  determine 
the  zeropoint  offsets  between  the  SDSS  and  Stetson’s  photometry  for  each  re- 
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u-g 

Figure  3.  An  illustration  of  the  effects  of  metallicity  on  photometric  trans¬ 
formations.  Here  the  u  —  g  color  serves  as  a  proxy  for  metallicity  (see  the 
right  panel  in  Fig.  2).  The  panels  shows  the  dependence  of  the  B  —  g  resid¬ 
uals  (top)  and  the  V  —  g  residuals  (bottom)  on  the  u  —  g  color.  The  effects 
of  metallicity  on  photometric  transformations  from  the  SDSS  to  the  BVRI 
system  appear  smaller  than  about  OP  01,  and  possibly  somewhat  larger  for  the 
B  —  g  transformation  for  stars  with  u  —  g  <  0.9. 


gion  separately  by  synthesizing  BVRI  magnitudes  from  SDSS  gri  photometry, 
and  comparing  them  to  Stetson’s  measurements.  The  implied  zeropoint  errors 
(which,  of  course,  can  be  due  to  either  SDSS  or  Stetson’s  dataset,  or  both)  are 
listed  in  Table  2.  For  regions  1-3  the  implied  errors  are  only  a  few  millimags 
(except  for  the  B  —  g  color  in  region  1 ) .  The  discrepancies  are  much  larger  for 
the  three  red  colors  in  region  4.  A  comparison  with  the  results  of  internal  SDSS 
tests  described  by  Ivezic  et  al.  (2006)  suggests  that  these  discrepancies  are  more 
likely  due  to  zeropoint  offsets  in  Stetson’s  photometry  for  this  particular  region, 
than  to  problems  with  SDSS  photometry.  We  contacted  P.  Stetson  who  con¬ 
firmed  that  his  observing  logs  were  consistent  with  this  conclusion.  Only  a  small 
fraction  of  stars  from  Stetson’s  list  are  found  in  this  region. 

Given  the  results  presented  in  this  Section,  we  conclude11  that  the  rms  for 
the  spatial  variation  of  zeropoints  in  the  SDSS  Stripe  82  catalog  is  below  01n01 
in  the  gri  bands. 


4.  Discussion  and  Conclusions 

Using  repeated  SDSS  measurements,  we  have  constructed  a  catalog  of  about 
700,000  candidate  standard  stars.  Several  independent  tests  suggest  that  both 
random  photometric  errors  and  internal  systematic  errors  in  photometric  zero- 
points  are  below  OP 01  (about  2-3  times  as  good  as  individual  SDSS  runs)  for 
stars  brighter  than  (19.5,  20.5,  20.5,  20,  18.5)  in  ugriz ,  respectively.  This  is  by 


11Here  we  assumed  that  it  is  a  priori  unlikely  that  the  SDSS  and  Stetson’s  zeropoint  errors  are 
spatially  correlated. 
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far  the  largest  existing  catalog  with  multi-band  optical  photometry  accurate  to 
~1%.  The  catalog  is  publicly  available  from  the  SDSS  Web  Site. 

In  this  contribution,  we  have  tested  the  photometric  quality  of  this  catalog 
by  comparing  it  to  Stetson’s  standard  stars.  Using  third-order  polynomials  to 
describe  the  color  terms  between  the  SDSS  and  BVRI  systems,  we  find  that 
photometric  measurements  for  main-sequence  stars  can  be  transformed  between 
the  two  systems  with  systematic  errors  smaller  than  a  few  millimagnitudes.  The 
spatial  variation  of  photometric  zeropoints  in  the  two  catalogs  typically  does  not 
exceed  0U01.  Consequently,  the  SDSS  Standard  Star  Catalog  for  Stripe  82  can 
be  used  to  calibrate  new  data  in  both  the  SDSS  ugriz  and  the  BVRI  systems 
with  a  similar  accuracy. 

Acknowledgments.  Funding  for  the  SDSS  and  SDSS-II  has  been  pro¬ 
vided  by  the  Alfred  P.  Sloan  Foundation,  the  Participating  Institutions,  the  Na¬ 
tional  Science  Foundation,  the  U.S.  Department  of  Energy,  the  National  Aero¬ 
nautics  and  Space  Administration,  the  Japanese  Monbukagakusho,  the  Max 
Planck  Society,  and  the  Higher  Education  Funding  Council  for  England.  The 
SDSS  Web  Site  is  http://www.sdss.org/.  The  SDSS  is  managed  by  the  As- 
trophysical  Research  Consortium  for  the  Participating  Institutions.  The  Partic¬ 
ipating  Institutions  are  the  American  Museum  of  Natural  History,  Astrophys- 
ical  Institute  Potsdam,  University  of  Basel,  Cambridge  University,  Case  West¬ 
ern  Reserve  University,  University  of  Chicago,  Drexel  University,  Fermilab,  the 
Institute  for  Advanced  Study,  the  Japan  Participation  Group,  Johns  Hopkins 
University,  the  Joint  Institute  for  Nuclear  Astrophysics,  the  Kavli  Institute  for 
Particle  Astrophysics  and  Cosmology,  the  Korean  Scientist  Group,  the  Chinese 
Academy  of  Sciences  (LAMOST),  Los  Alamos  National  Laboratory,  the  Max- 
Planck-Institute  for  Astronomy  (MPIA),  the  Max-Planck-Institute  for  Astro¬ 
physics  (MPA),  New  Mexico  State  University,  Ohio  State  University,  University 
of  Pittsburgh,  University  of  Portsmouth,  Princeton  University,  the  United  States 
Naval  Observatory,  and  the  University  of  Washington. 


References 

Abazajian,  K.,  Adelman,  J.  K.,  Agueros,  M.,  et  al.  2003,  AJ,  126,  2081 
Abazajian,  K.,  Adelman,  J.  K.,  Agueros,  M.,  et  al.  2004,  AJ,  128,  502 
Abazajian,  K.,  Adelman-McCarthy,  J.  K.,  Agueros,  M.,  et  al.  2005,  AJ,  129,  1755 
Adelman-McCarthy,  J.  K.,  et  al.  2006,  ApJS,  162,  38 
Allende  Prieto,  C.,  Beers,  T.  C.,  Wilhelm,  R.,  et  al.  2006,  ApJ,  636,  804 
Fukugita,  M.,  Ichikawa,  T.,  Gunn,  J.  E.,  et  al.  1996,  AJ,  111,  1748 
Gunn,  J.  E.,  Carr,  M.,  Rockosi,  C.,  et  al.  1998,  AJ,  116,  3040 
Gunn,  J.  E.,  Siegmund,  W.  A.,  Mannery,  E.  J.,  et  al.  2006,  AJ,  131,  2332 
Hogg,  D.  W.,  Finkbeiner,  D.  P.,  Sclilegel,  D.  ,J.  &  Gunn,  J.  E.  2002,  AJ,  122,  2129 
Ivezic,  Z.,  Lupton,  R.  H.,  Anderson,  S.,  et  al.  2003,  in  Variability  with  Wide  Field 
Imagers,  Mem.  Soc.  Ast.  It.,  74,  978  (also  astro-ph/0301400) 

Ivezic,  Z.,  Lupton,  R.  H.,  Schlegel,  D.  J.,  et  al.  2004,  AN,  325,  583 
Ivezic,  Z.,  Smith,  J.  A.,  Miknaitis,  G.,  et  al.  2006,  submitted  to  AJ 
Landolt,  A.U.  1983,  AJ,  88,  439 
Landolt,  A.U.  1992,  AJ,  104,  340 

Lupton,  R.  H.,  Gunn,  J.  E.,  &  Szalay,  A.  1999,  AJ,  118,  1406 

Lupton,  R.H.,  Gunn,  J.  E.,  Ivezic,  Z.,  et  al.  2001,  in  ASP  Conf.  Ser.  238,  Astronomical 
Data  Analysis  Software  and  Systems  X,  eds.  F.  R.  Harnden,  Jr.,  F.  A.  Primini, 


Comparison  of  SDSS  Standards  with  Stetson ’s  Standards 


175 


and  H.  E.  Payne  (San  Francisco:  Astr.  Soc.  Pac.),  p.  269  (also  astro-ph/0101420) 
Lupton,  R.  H.,  Ivezic,  Z.,  Gunn,  J.  E.,  et  al.  2003,  in  “Survey  and  Other  Telescope 
Technologies  and  Discoveries”,  Tyson,  J.  A.  &  Wolff,  S.,  eds.  Proceedings  of  the 
SPIE,  4836,  350 

Oke,  J.  B.,  &  Gunn,  J.  E.  1983,  ApJ,  266,  713 
Pier,  J.  R.,  Munn,  J.  A.,  Hindsley,  R.  B.,  et  al.  2003,  AJ,  125,  1559 
Schlegel,  D.,  Finkbeiner,  D.  P.  &  Davis,  M.  1998,  ApJ500,  525 
Scranton,  R„,  Johnston,  D.,  Dodelson,  S.,  et  al.  2002,  ApJ,  579,  48 
Sesar,  B.,  Svilkovic,  D.,  Ivezic,  Z.,  et  al.  2006,  AJ,  131,  2801 
Smith,  J.  A.,  Tuckey  D.  L.,  Kent,  S.  M.,  et  al.  2002,  AJ,  123,  2121 
Smolcic,  V.,  Ivezic,  Z.,  Gacesa,  M.,  et  al.  2006,  MNRAS,  371,  121 
Stetson,  P.  B.  2000,  PASP,  112,  925 
Stetson,  P.  B.  2005,  PASP,  117,  563 

Stoughton,  C.,  Lupton,  R.  H.,  Bernardi,  M.,  et  al.  2002,  AJ,  123,  485 
York,  D.  G.,  Adelman,  J.,  Anderson,  S.,  et  al.  2000,  AJ,  120,  1579 


Discussion 

Stubbs:  Can  you  use  the  high-extinction  data  (taken  through  clouds)  to  extend 
the  bright  limit  to  the  TYCHO  catalog? 

Ivezic:  One  needs  about  3  magnitudes  of  extinction  to  extend  the  SDSS  satura¬ 
tion  limit  to  the  TYCHO  faint  limit.  We  have  some  data  with  cloud  extinction 
larger  than  3  mags,  but  not  much.  It  may  be  enough  to  find  a  few  stars  and  to 
directly  compare  the  two  photometric  systems.  Thus,  we  need  more  bad  data! 

Landolt:  If  one  sees  a  point  of  interest  in  one  of  your  color-color  plots,  can  you 
identify  the  star  itself? 

Ivezic:  Yes,  one  can  search  the  list  by  provided  color  coordinates,  and  then  look 
up  the  sky  coordinates. 

Henden:  Your  plots  indicate  good  variable-constant  source  separation.  What  is 
the  fraction  of  variable  sources  that  you  find? 

Ivezic:  Depends  on  the  cutoff.  For  5%  rms  cutoff  in  the  g  and  r  bands,  about 
10%  of  point  sources  are  variable.  However,  the  precise  value  is  very  sensitive 
to  the  adopted  cutoff  in  magnitude.  A  fair  statement  would  be  “between  5  and 
20%  of  faint  optical  point  sources  are  variable  at  the  5%  rms  level” . 

Ageorges:  How  good  is  the  astrometric  precision  of  the  SDSS  catalog? 

Ivezic:  The  astrometric  accuracy  of  single  SDSS  scans  at  the  bright  end  is  0,/l 
(absolute)  and  20-30  mas  (relative).  The  “multi-epoch”  catalog  is  expected 
to  be  better,  similarly  to  photometry,  but  we  don’t  know  yet  by  how  much. 
The  expectation  is  about  10  mas  relative  accuracy  per  coordinate  for  the  mean 
position. 


